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Abstract. The purpose of this paper is to show how the problem of finding the zeros 
of unilateral n-order quaternionic polynomials can be solved by determining the eigen- 
vectors of the corresponding companion matrix. This approach, probably superfluous in 
the case of quadratic equations for which a closed formula can be given, becomes truly 
useful for (unilateral) n-order polynomials. To understand the strehgth of this method, 
we compare it with the Niven algorithm and show where this (full) matrix approach 
improves previous methods based on the use of the Niven algorithm. For the convenience 
of the readers, we explicitly solve some examples of second and third order unilateral 
quaternionic polynomials. The leading idea of the practical solution method proposed 
in this work can be summarized in following three steps: translating the quaternionic 
polynomial in the eigenvalue problem for its companion matrix, finding its eigenvectors, 
and, finally, giving the quaternionic solution of the unilateral polynomial in terms of the 
components of such eigenvectors. A brief discussion on bilateral quaternionic quadratic 
equations is also presented. 

Subject Classification: 02.10.Yn - 02.10.Yn (PACS), 15A33 - 20G20 (MSG). 



I. INTRODUCTION 



As far as we knovif, the problem of solving quaternionic quadratic equations and the study of the 
fundamental theorem of algebra for quaternions were first approached by Niven and Eilenbcrg [1,2]. 
After these fundamental works, the question of counting the number of zeros of n-order quaternionic 
polynomials and the problem to find the solutions have been investigated. An interesting application 
of unilateral quadratic equations (where a closed formula for their solutions can be obtained) is 
found in solving homogeneous quaternionic linear second order differential equations with constant 
coefficients [3,4]. The solution of 

— ^{x)~a,—-^{x)-a,^{x)^Q, (1) 
da;'^ da; 

where 4* : K H , ao.i G H (the set of quaternions) and a; e M, can indeed be reduced, by setting 
'^{x) = exp[(7a;] {q E H) and using the H-linearity of Eq.(IT]), to the following quadratic equation 



— cii q — ciq = 



(2) 
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Generalizing this result, we can see that the problem of finding the solution of ?i-ordcr (H-linear) 
differential equations with constant coefficients, 

d" d*" 

where ao.i....,n-i ^ HI, can be immediately transformed to the problem of finding the zeros of the 
corresponding n-order unilateral quaternionic polynomial, 

n-i 

A„(g):=g"-Ea.g^ . (4) 

s=o 

The important point to be noted here is that such a solution method based on finding the zeros of the 
corresponding quaternionic polynomial for H-linear differential equations with constant coefficients 
does not work for M and C-linear differential equations. For example, in non-relativistic quaternionic 
quantum mechanics [5] , the one-dimensional Schrodinger equation (for quaternionic stationary states) 
in presence of a constant quaternionic potential, reads 

i — —^x)-{iV,+jV, + kVs) <P{x)^^Mx)Ei , (5) 
Zm ax 

and, due to the C-linearity, its solution cannot be expressed in terms of a quaternionic exponential 
function, exp[ga;]. Recent studies of quaternionic barrier [6] and well [7] potentials confirmed that the 
solution of Eq.([5|) has to be expressed as the product of two factors: a quaternionic constant coefficient, 
p, and a complex exponential function, exp[zx] (z € C). Using $(2;) = p exp[zx], the Eq.Q reduces 
to 

i—pz^ - {iV,+ jV2 + kV,) p^ -pEi , (6) 

which obviously does not represent a quaternionic unilateral polynomial. For a more complete dis- 
cussion of M and C linear quaternionic differential equation theory, we refer the reader to the papers 
cited in rcfs. [3,8]. 

The problem of finding the zeros of unilateral quaternionic polynomials was solved in 1941 by 
using a two-step algorithm. In his seminal work [9], Niven proposed to divide the unilateral n-order 
quaternion polynomial jl]) by a quadratic polynomial with real coefficients (co,i S M) 

C2{q) ■.= q- -CiQ-Co . (7) 

Then, by using the polynomial equation 

Mq)=B,-,_,iq)C,iq)-D,{q) , (8) 

where 

n — 3 

Bn-2{q) g""' - ^ b, q' , D,{q) --d.q + do , bo,i,...n-3 , do,t G H , 

the problem of finding the zeros of A„ (g) was translated in the following two-step problem: 

— step 1 - determine and d^ in terms of Co,i and ao_i_..._„_i; 

— step 2 - obtain two coupled real equations which allow to calculate Cq and Ci. 

To facilitate the understanding of this paper, we explicitly discuss the Niven algorithm in the next 
section. Since its proof is quite detailed and, furthermore, since its use can be completely avoided, we 
merely state the results without proof. 

In a recent paper [10], Serodio, Pcrcira and Vitoria (SPV) proposed a practical technique to 
compute the zeros of unilateral quaternion polynomials. The SPV approach is based on the use of the 
Niven algorithm [9]. Their idea was to improve such an algorithm by calculating the real coefficients 
Co and Ci by using, instead of the second part of the Niven algorithm, the (complex) eigenvalues of the 
companion matrix associated with the quaternionic polynomial An{q). The SPV approach, avoiding 
the solution of the coupled real equations (step 2), surely simplifies the Niven algorithm and gives 
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a more practical method to find the zeros of n-order unilateral quaternionic polynomials. However, 
once the real coefficients Cq and Ci are obtained in terms of the moduli and real parts of the complex 
eigenvalues of the companion matrix, we have to find the quaternionic coefficients do and di . It seems 
that we have no alternative choice and we have to come back to the Niven algorithm (step 1). 

In this paper, we show that the matrix approach (based on the use of the quaternion eigenvalue 
problem) leads us to find directly the solutions of unilateral quaternion polynomials by calculating 
the eigenvectors of the companion matrix. This allow us to completely avoid the use of the Niven 
algorithm and, consequently, to hardly simplify the method to find the zeros of unilateral quaternionic 
polynomials. 

We interrupt our discussion at this point to introduce, the Niven algorithm and the vector-matrix 
notation. These tools would permit a more pleasant reading of the article and rend our exposition 
self-consistent. 



II. THE NIVEN ALGORITHM 

In this section, we aim to concentrate on the Niven algorithm. We shall first discuss the method 
proposed by Niven for a generic n-order quaternionic polynomial, then, to illustrate the potential and 
problems of the Niven algorithm, we consider an explicit example, i.e. a quadratic polynomial. 

In order to determine the quaternionic coefficients dg and c?i in terms of the real coefficients of the 
second order polynomial C2{q) and of the quaternionic coefficients of the n-order polynomials An{q), 
we expand the polynomial product _B„_2 ((7)6*2 (g), 

5n-2(g)C2(9) = g" - (6„_3 + ci) g"-' - (5„_4 - c, 5„_3 + Co) g""' + 

n-3 

- ^ (6s-2 - Ci - Cn bs) q" + 

S = 2 

+ (ci ba + Co bi)q + Co 5o . (9) 
For a particular choice of the quaternionic coefficients of the polynomial Bn--2{q)^ i.e. 

bn-3 — Q-n-l Ci , 

bn-4 ~ 4" Ci bji^^ Co , 

bs_2 ^ as + Cibs^i + Cobs (s = 2, 3, • • • , n - 3) , (10) 
we can rewrite Eq.® as follows 

-B„_2 (9)6*2(9) = An{q) + (fli -f Ci 60 + Co 61) g -I- ao -I- Co 60 ■ (H) 
Comparing this equation with the polynomial equation ([8|), we obtain (step 1) 

di = CLi -\- Ci feo(Co,i; (i2,3,--- ,n-l) ^~ Co ^i(Co,i5 (i2,3,--- ,n-l) : 

do = flo + Co &o(Co,i; a2,3,--- ,n-i) ■ (12) 

This show that by some simple algebraic manipulations, we can express the quaternionic coefficients 
do,i in terms of the real coefficients Co,i and of the quaternionic coefficients ao,i,...,n_i. This represents 
step 1 of the Niven algorithm. As stated in the introduction, to complete the Niven algorithm, we have 
to determine the real coefficients Co.i (step 2). Once such coefficients arc obtained, we can explicitly 
calculate the coefficients do^i by using Ea. (|12p which, after step 2, only contains known quantities. It 
is natural to ask why the quaternionic coefficients c?o,i are important in calculating the quaternionic 
solution of the polynomial An{q). The answer is given by observing that if q^, is a solution of An{q) 
An{q^) = 0] and { Co , Ci } = {- , 2Re[g*] } 6*2(9*) = 0], then from Eq.®, we find 

A(g*) = => q, = -d,do/\d,\' . (13) 

The quaternionic solution can thus be expressed in terms of the coefficients do,i. Turning to the 
problem of determining the real coefficients Co,i, Niven noted that 



Co = -k*r = -Mor/Mir and ci= 2Rc[g,]--2Rc[Jido]/Mir 
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Consequently, in order to find Co,i we have to solve the following system (step 2) 

Co I'^l (Co,i i <^l,2,3,--- ,n-l) I ~t" |<io(*^o,li ^0,2,3, ■■■ ,n-l)| =0 , 
Ci |di(Coa;ai,2,3,...,„_i)|^ + 2 Rc[Ji(Co,i; 01,2,3, •■■ ,n-i)do(Co,i; Ao, 2,3, •■■ ,n-i)] = . (14) 

Each real coupled solution { Co , Ci } gives the modulus and the real part of the quaternionic polyno- 
mial solution. 

• Second order quaternionic polynomials. The main (practical) problem in using the Niven 
algorithm is in step 2, that is finding the real solutions of the coupled system (fT4)) . To illustrate the 
difficulty in using the Niven algorithm, we solve the following quaternionic quadratic equation 

+ j q + - k) = [^2(9) = - Oi g - Oo , flo = fc - 1 , ai = -.?]. (15) 

Second order polynomials represent the more simple situation in which we can test the Niven algo- 
rithm. In this case, the quaternionic polynomial B^-^iq) reduces to B„{q) = —b„ ~ 1. Consequently, 

Ba{q)C2{q) = - Ci q - Co 

= - fli q - flo + ( fli - Ci ) q -f flo - Co 
= A,iq)+D,iq) . 

From the previous equation or equivalently by setting = and 5o = — 1 in Eq. ()12p . we obtain 

rfi = fli - Ci , 

do = Co - Co . (16) 

Now, for second order polynomials. Eg. pi]) reduces to 

Co |ai - Ci|^ + \ao - Col'' = i 
c, \a, - c,\' + 2Re[(ai - cO(a„ - c„)] = . (17) 

The solution of this system is not simple. In this special case, the choice of 

tti = —j and tta ~ k — I , 

reduces the system ([T7|) to 

Co (1 + ci) + (1 + c,y + 1 = 0, 

Ci(l + cJ) + 2ci(l + Co) = . (18) 
The discussion can be then simplified by considering 

Ci 7^ and Ci — . 
In the first case no real coupled solution exists. For Ci = 0, we have 

(c c )-/(-i'0) ^ ^ ^ r(fc,-j) 



Finally, by using Eq. p3|) we conclude that the solutions of the quadratic equation p5|) are given by 

By direct substitution, we can easily verify that —i and —{i + j) represent the zeros of our quadratic 
polynomial. Before concluding this brief review of the Niven algorithm and underlining that it surely is 
an important (old) technique for obtaining the solution of n-order unilateral quaternionic polynomial, 
we wish to emphasize that the system (|17p and, in general, the system (|12|) are not very practical and 
their solutions could require laborious calculations. 
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III. COMPANION MATRIX AND SPV METHOD 



Taking as our guide the standard complex matrix theory, we associate to the n-ordcr unilateral 
quaternionic polynomial An{q) the following companion matrix 



1 
1 



V 







•• 





•• 





•• 1 


0/ 



(20) 



In the standard complex theory, the zeros of n-order polynomials can be determined by calculating the 
eigenvalues of the corresponding companion matrix. Any direct attempt based upon a straightforward 
expansion of the determinant |Af„[H] — A/| is surely destined for failure due to the non-commutativity 
of the quaternionic algebra [11]. The quaternionic eigenvalue problem has recently been investigated 
by means of the complex representation of quaternionic matrices. It is clear that the rigorous presen- 
tation of the quaternionic eigenvalue problem would require a deep knowledge of quaternionic matrix 
theory. To facilitate the understanding of this paper, we do not discuss any theoretical aspect of 
the quaternionic eigenvalue problem here and aim to give a practical matrix solution method based 
on complex translation rules. We refer the reader to the articles cited in refs. [12-14] for a through 
discussion of quaternionic diagonalization and Jordan form. In order to appreciate some of the poten- 
tialities of the complex representation of quaternionic matrices, we shall consider in detail the case 
of second-order polynomials and show how we can immediately obtain the coefficients £„ and Ci by 
calculating the real part and the modulus of the eigenvalues associated to the complex matrix which 
represents the complex translation of the quaternionic companion matrix (SPV method [10]). 

Before analyzing the companion matrix associated to a second-order unilateral quaternionic poly- 
nomial and calculate its eigenvalues, we give the main idea of the complex translation by using a 
generic 2x2 quaternionic matrix. 



AIM] = 



9i 92 

93 94 



9l,2,3,4 € 



The complex translation would proceed along the following lines. A first possibility is the use of the 
symplectic form directly for the quaternionic matrix A'f2[E[], i.e. 



M,[H] 





Zi \ 








+J 






Zi ) 






(21) 



where ^1^2,3,4 and 1^1,2,3,4 G C. In this case the symplectic decomposition directly applies to the complex 
matrix blocks. Another possibility is represented by the use of the symplectic decomposition for the 
quaternionic elements of M2[IHI], i.e. 



M2[H] 



Zi+ JWi Z2+ J W2 
23 + j W3 Z4 -I- j Wi 



M,[C\ 




(22) 



For the problem of finding the zeros of quaternionic polynomial, the complex representations (|21[) and 
(|22p represent equivalent choices. In fact, there exist a permutation matrix 



P. 



/ 1 \ 

10 

10 

V 1 y 



which convert M4[C] to M4[C] by using a simple similarity transformation, 

MAC] ^p^MAap^ . 
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In Appendix A, we extend this similarity for the matrices A/2n[C] and Af2n[C] (complex counterparts 
of the n-dimension quaternionic matrix AIn [H] ) , 

by explicitly giving the permutation matrix Paji • Such a matrix will be obtained by working with the 
complex representations of the quaternionic column vectors associated to the matrix representations 
(|2ip and (|22p . Due to this equivalence, in all the explicit examples given in this paper, we shall use, 
without loss of generality, the complex matrix representation A/an [C] . 

• Second order quaternionic polynomials. Let us consider the unilateral quaternionic polynomial 
discussed in the previous section, Eq. (|15[) . In this case, the quaternionic companion matrix is 



AUm-{^ ^ \ ^ ) • (23) 

The SPV method allow to obtain the real coefficient Cf, and Cj by calculating the eigenvalues of the 
complex matrix M4 [C] , 

{ Ai , A2 , Ai , A2 } , 

which, due to the particular structure of the complex symplectic translation, always appear in con- 
jugate pairs [12]. By using the complex translation rules (|21[) . we immediately obtain the following 
complex counterpart of A/2P], 

0-1 1-i 

w = I _ ■ h (24) 
0010 

whose eigenvalues can be easily calculated by using jAfjp] — A/4I = 0. As observed above, the 
eigenvalues appear in conjugate pairs, 

{ Ai = i , A2 = i \/2 , Ai = - i , A2 = - I \/2 } . 

The real coefficients Cq and Ci, 



( Co , Ci ) 



(-|Aih2Re[Ai]) = (-l,0) 
(-|A2h2Re[A2]) = (-2,0) 



can now be obtained without the need to solve the system ^TE\i . This method proposed by Serodio, 
Pereira and Vitoria [10] avoids the solution of the complicated coupled system ([T4| and represents a 
simpler way to find the real coefficients Cq and Ci. Consequently, it improves the Niven algorithm. Note 
that step 1 in the Niven algorithm has yet to be used to find the quaternionic coefficients do Sund di . 
In the next section, we shall prove that by calculating the eigenvectors of the complex representation 
of the quaternion companion matrix Af„[IHI], we can directly obtain the zeros of 71-order unilateral 
quaternionic polynomials. 



IV. FINDING THE ZEROS BY USING THE EIGENVECTORS METHOD 

In the complex theory, the zeros of the n-order polynomial, j4„(z) = 0, are determined by calculating 
the eigenvalues of the corresponding companion matrix, Af„[C], i.e. |A^„[C] — A/| =0. For quater- 
nions, this method cannot be carried through successfully. In fact, the eigenvalues of the quaternionic 
companion matrix are not univocally determined. Consequently, we have to pursue a different track. 
Before demonstrating two theorems which will guide us to the solution of our quaternionic math- 
ematical challenge, wc briefly recall the main results concerning the (right) eigenvalue problem for 
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quaternionic matrices. Due to the non-commutativity of the quaternionic algebra, the eigenvalues of 



M„[H] 



1 • 
1 • 



V 












/ 



/ {j,k) {j.k) {j.k) {j.k) \ 



V 



/ 



(25) 



W„[C] 



belong to the following class of equivalence 

where (by using the standard definition adopted in literature) A„ are the eigenvalues of 



M,„[C] = 



W„[C] 



z*\ 



(26) 



(27) 



with Iin[A„] > 0, and u„ are quaternionic unitary similarity transformations. The starting point in 
finding the zeros of n-order quaternionic unilateral polynomials is the same as in the complex theory, 
i.e. the solution of the eigenvalue problem for the companion matrix. 



M„[M] $ = $ A , A e C , Im[A] > , 



(28) 



and 



^2 



Via.,....n e 



Nevertheless, once the complex eigenvalues have been calculated we also get infinite similar quater- 
nionic eigenvalues 

A/„[M] = fwAw , (29) 

where 4' = $ u, u G H, and uu = 1. Geometrically, the quaternionic eigenvalue uXu represents, for 
the imaginary part of A, the following three-dimensional rotation: 

A 

i 

(Im[A],0,0) 

/\ 

6* = 2 arctan[|(Re[— iu], Re[— ju], Re[— fcu])|/ito] , axis: (Re[— m], Re[— j'w], Rc[— fcu]) 

\ / 

(Re[— mAw], Re[— juAw], Re[— fcuAu]) 

i 

uXu . 

We demonstrate later that, for each complex eigenvalue. A, it is possible to choose an appropriate 
similarity transformation u*-*^ Au*-*^ which corresponds to one of the zeros of the n-order quaternionic 
unilateral polynomial A„(q). This breaks down the symmetry between the infinitely many equivalent 
directions in the eigenvalue quaternionic space, uXu. The complete set of the zeros of the n-order 
quaternionic unilateral polynomial An{q) will be then given by 



(30) 



Now, we proceed with the main objective of this Section and prove that the eigenvalue spectrum ([30|) 
is connected with the zeros of An{q). In doing that, the first step will be the proof of two theorems. 
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Theorem A and Theorem B concern the (right) complex eigenvalue problem (j28|) . 

Theorem A. The last component of the quaternionic eigenvector $ is non null. 

Proof. Due to the special form of the companion matrix, from Eq. (pS)) wc get the following equalities: 

fln-l fl + a„-2 1- fll fn-l + flo <<5n = <y5i A , 



(31) 



Thus, (p„ = implies </5i,2,...,n-i — 0. This means 

(y5„ = =^ $ = ■ 

Theorem B. Let <I> he an eigenvector of the companion matrix M„[]HI] associated with the polynomial 
An{q). The solution q, expressed in terms of the two last components of^, is given by fn-ifn^- 

Proof. Combining the last n — 1 equalities in (|3ip . we obtain 

"^n-l = <Pti A , 

'Pn-2 = <Pn-l A = A^ , 



This implies 



\'~Pn/ 



<^n A" 



/A"-\ 



V 1 / 



The first equation of system (PT|) can be then re-written as follows 

fln-i "ySn A"~^ + a„_2 "y^n A"~^ + • • • + flj (y9„ A + Cq iy9„ = iy9„ A" . 

Theorem A guarantees that ipn ^ 0. We can thus multiply (from the right) the previous equation by 
(/'^^ obtaining 



fln-l A" Ifn + '^n~2 <y5„ A" H h flj A + Oq = (y5„ A" i^. 



n 7 



or equivalently 

a„_i ( (y9„ A )""' + a„_2 ( A (p;;' )"~^ H h fli A (^s;;' + Oo = ( (p„ A (^s,";' )" . 

Consequently, g = A(p„^ and by using that ipn-i = VnA 

g = v3„-i(p„' . (32) 

This completes the proof. ■ 

The results of Theorem A and Theorem B can be easily extended to the (right) quaternionic eigenvalue 
problem (|^^ . hi fact. 



and 



q^ipn-ifn => q = iPn-lU{ipnU) (/3„_i U M (^„' = ■071-1 V'n ' 



S. De Leo, G. Ducati, V. Leonardi: unilateral quaternionic polynomials 



9 



We conclude this subsection going back to the discussion on the quaternionic similarity transformation 
which breaks down the symmetry between the infinitely many equivalent directions in the eigenvalue 
quaternionic space. By using Theorem B, wc can write the zero of the quaternionic polynomial An{q) 
in terms of the last component of the eigenvector $ and of the complex eigenvalue A, 

Q = fnXfn' ■ (33) 

This allows us to define (for the companion matrix) a privileged (right) quaternionic eigenvalue prob- 
lem, i.e. 

M„[H] = u(«) A (u'«) = ipn/\iPn\ ) 

= q , (34) 



with 



^(9) ^^if,-' |(^„| 



V 1 / 



Recall that 



ipn-i^Pn^ and A = <^„V«-i- 



For the complex case, commutativity guarantees q = X. The problem of finding the zeros of n- 
order unilateral quaternionic polynomials is enormously simplified by using the results obtained in 
this Section. This will be explicitly seen in the next subsections, where the solutions of second and 
third order quaternionic unilateral equations are calculated using complex matrix translation. For the 
convenience of the readers, we illustrate the steps of our method in the following sequence: 

An{q) (n-order unilateral quaternionic polynomial) , 

Af„[]HI] (quaternionic companion matrix) , 

M2„[C] (complex translated matrix) , 

(wi, Wa, uJm cTi, (72, <7n)\ (eigenvector of M2„[C] corresponding to the complex eigenvalue A) 

{lui + j ai, uj^ + j uin + j a'n)\ (quaternionic translated eigenvector $>) , 

q ~ ifin-i Vn^ (zero corresponding to A) , 

A ~ ip^^ <fn-i (check of the solution) . 

The importance of the new method is due, in great extent, to the possibility of using complex matrix 
analysis. 



IV-A. SECOND ORDER UNILATERAL POLYNOMIALS 



In Section II and in Section III we have solved (by means of the SPV algorithms) the following 
(unilateral) quaternionic quadratic equation 



g'+j9 + (i-fc) = 



10 
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We now calculate the zeros of A2 (q) by using our method: 



A,{q)=q'+jq+il-k) , 




LJ2 + j (J2 



1 



x/2 . 



/V2(V2 - 1 - fc) 
j - *(V2-1) 



;V2 . 



IV-B. THIRD ORDER UNILATERAL POLYNOMIALS 



Let us now give an example of solution for a third order quaternionic polynomial. 
^3(9) ^q^ + kq^ +iq-j , 

MM] = 




Me[C] 



/ 0-i i 

1 

10 

1 1 i 
1 

yo 1 



-1\ 






0/ 




V2 



i + 1 

i - 1 
\ V2 J ^ 



- 1 - kV2\ 
V2 + j-k 

k~v2 \ r , 



-V2-j-k 
l + t+jV2 J 



V2 
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V. BILATERAL QUADRATIC EQUATIONS 

So far, we have not analyzed the possibihty to have left and right acting coefficients in the quaternionic 
polynomials. This topic exceeds the scope of this paper. At the moment it represents an open ques- 
tion and satisfactory resolution methods have to be further investigated. Nevertheless, an interesting 
generalization of the matrix approach introduced in the previous section can be found in discussing 
bilateral equations [15], i.e. 

p'^ — ai p + p Pi - ao ~ . (35) 

For second order bilateral quaternionic polynomial, we introduced the following generalized companion 
matrix 

MM^I^'l' 2) ' (36) 
which reduces to the standard one for f3i ~ 0. The (right) complex eigenvalue problem for this matrix, 

7W2[IHI]$ = $A, (37) 

is equivalent to 

where (pi^2 arc the components of the eigenvector (f>. It is immediate to show that = implies 
$ = 0. Consequently, the second component of the (non trivial) eigenvector $ is always different from 
zero. Multiplying the first of the previous equations by (^,7^ (from the right), we obtain 

ai ipi ip^^ + cko = if^' ip2 A (fi^' . 

From the second equation, we get 
Finally, 

2 

which, once compared with Eq. (j35p . implies 

P = Vi ^2^ ■ (38) 

The solution for the bilateral quadratic equation (j35p is then formally equal to the solution given for 
the unilateral one. The complex eigenvalue A can be expressed in terms of the $ components and of 
the right acting coefScient /3i by 

A = <^JVi + </'^'A'P2 • (39) 



It is worth pointing that the bilateral equation (j35[) can be reduced to the following equivalent uni- 
lateral equation 

- fli q - ao = , (40) 
where /?! by using 

p^q~Pi. (41) 

The next examples should emphasize one more time the advantage of using of the matrix approach 
proposed in this paper to solve unilateral (and some particular bilateral) quaternionic equations. It 
would be desirable to extend this approach to more general polynomials (with left and right acting 
quaternionic coefhcients) but we have not been able to do this and from our point of view this 
should deserve more attention and further investigation. The question of solving general quaternionic 
polynomials is at present far from being solved. 
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V-A. SECOND ORDER BILATERAL POLYNOMIALS 

Let us consider the bilateral quadratic equation 

p^ — ip + pj — k = 



(42) 



To solve this equation, we shall follow the same steps of the resolution method presented in the 
previous section: 

p^ ~ip + pj - k , 
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V-B. EQUIVALENT SECOND ORDER UNILATERAL POLYNOMIALS 

In this subsection, we solve the unilateral quadratic equation obtained from Eq. (|42p by using p = 9— j: 



i - 1 

10 

1 -i 

10 




/-^/2(l + ^/2)\ 
i (1 + \/2) 

1 



-x/2(l + V2 + fc)\ 
z(l + V2)+j 



g = <^,^-i = {0, i + j} A = ^->i ={0, i\/2}] 

As expected, by recalling that 

P = 9-/3i=g-j = {-J,i} 
we recover the solution of Eq. obtained in the previous subsection. 
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VI. CONCLUSIONS 

In this paper we have presented a practical method to find the roots of unilateral quaternion polyno- 
mials by means of the eigenvalue spectra and eigenvectors of the companion matrix. In particular, we 
have shown that the zeros of unilateral quaternion polynomials can be directly given in terms of the 
two last component of the eigenvectors of the companion matrix. The practical resolution method is 
enormously simplified by using the complex matrix translation. We have also extended the method 
to bilateral equations introducing the generalized companion matrix. One question still unanswered 
is whether similar approaches could be used to finding the zeros of n-ordcr polynomial with left and 
right acting quaternionic coefficients. The paper is essentially self-contained and the possibility to 
translate the quaternionic eigenvalue problem to the corresponding complex one, then, come back to 
quaternions, and give the final solutions should represent a fundamental tool for the reader which 
is not so familiar with the quaternionic algebra. In this spirit, the applications (explicit examples 
of 2 and 3 order quaternionic polynomials) would clarify the advantage in adopting this full matrix 
approach in solving quaternionic unilateral equations. 

ACKNOWLEDGEMENTS. The authors wish to express their gratitude to Prof. Nir Cohen for 
several helpful comments and for many stimulating conversations. One of the author (S.D.L) greatly 
thanks the Department of Mathematics (Parana University), where the paper was prepared, for the 
invitation and hospitality. This is part of the research project financed by the Araucaria Foundation 
(V.L and G.D). 



References 

1. I. Niven, "The roots of a quaternion" , Amer. Math. Monthly 49, 386-388 (1942). 

2. S. Eilenberg and I. Niven, "The fundamental theorem of algebra for quaternions". Bull. Amer. Math. Soc. 
50, 246-248 (1944). 

3. S. De Leo and G. Ducati, "Quaternionic differential operators", J. Math. Phys. 42, 2236-2265 (2001). 

4. S. De Leo and G. Ducati, "Solving simple quaternionic differential equation", J. Math. Phys. 43, 2224- 
2233 (2003). 

5. S. L. Adler, Quaternionic quantum mechanics and quantum fields, (New York: Oxford University Press, 
1995). 

6. S. De Leo, G. Ducati and C. C. Nishi, "Quaternionic potentials in non-elativistic quantum mechanics", 
J. Phys. A 35, 5411-5426 (2002). 

7. S. De Leo and G. Ducati, "Quaternionic bound states", J. Phys. A 38, 3443-3454 (2005). 

8. S. De Leo and G. Ducati, "Real linear quaternionic differential equations", Comp. Math. Apll. 48, 1893- 
1903 (2004). 

9. I. Niven, "Equations in quaternions", Amer. Math. Monthly 48, 654-661 (1941). 

10. R. Serodio, E. Pereira and J. Vitoria, "Computing the zeros of quaternion polynomials", Comp. Math. 
Appl. 42, 1229-1237 (2001). 

11. N.Cohen and S. De Leo, "The quaternionic determinant". El. J. Lin. Alg. 7, 100-111 (2000). 

12. S. De Leo and G. Scolarici, "Right eigenvalue equation in quaternionic quantum mechanics", J. Phys. A 
33, 2971-2995 (2000). 

13. S. De Leo, G. Scolarici and L. Solombrino, "Quaternionic eigenvalue problem", J. Math. Phys. 43, 5812- 
2995 (2002). 

14. T. Jiang, "Algebraic methods for diagonalization of a quaternion matrix in quaternion quantum theory", 
J. Math. Phys. 46, 052106-8 (2005). 

15. R. Carvalho and S. De Leo, "Quaternions in Algebra and Analysis" (Campinas: DMA-Unicamp, 2002). 



14 



S. De Leo, G. Ducati, V. Leonardi: unilateral quaternionic polynomials 



APPENDIX A. EQUIVALENCE BETWEEN TRANSLATION RULES 



Let 



+ j 



/.A 

(72 



(To 



(43) 



and 



(^2 + j 



0-1 
<T, 



(44) 



be the complex vectors obtained by using the complex translation rules given in Eqs. (|2ip and (p2l 
A simple algebraic calculation shows that 



CTi 
(J2 



(jJr 



— P2n 



UJ„ 



(45) 



where 



1 (r,s) = (1,1), (2,71+1), (3,2), (4,71 + 2), (2n-l,n), (2n,2n) , 

otherwise . ^ ' 



The foregoing result can be then used to prove the equivalence between the complex matrices Man [C] 
and il/2„[C], 

M2n[q =P2„Af2„[C]P2„ . (47) 



